A numerical approach to model chemistry of complex organic molecules in a protoplanetary disk

Abstract Multiphase astrochemical modeling presents a numerical challenge especially for the simulation of objects with wide range physical parameters such as protoplanetary disks. We demonstrate an implementation analytical Jacobian integration system differential rate equations that govern chemical evolution in star-forming regions. The allowed us to greatly improve stability code disk conditions. utilize MONACO study abundances species model includes 670 and 6,015 reactions gas phase on interstellar grains. specific feature utilized is inclusion low-temperature processes leading formation complex organic molecules (COMs), included previously models chemistry COMs prestellar clouds. To test impact simulations disks, we calculated composition using two-phase four variants reaction network, three values surface diffusion rates, two types initial composition. also show preliminary three-phase model.